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Abstract. We have obtained the full 1-200 fim spectrum 
of the low luminosity (36 Lq) Class I protostar Elias 29 
in the p Ophiuchi molecular cloud. It provides a unique 
opportunity to study the origin and evolution of interstel- 
lar ice and the interrelationship of interstellar ice and hot 
core gases around low mass protostars. We see abundant 
hot CO and H2O gas, as well as the absorption bands of 
CO, CO2, H2O and “6.85 /xm” ices. We compare the abun- 
dances and physical conditions of the gas and ices toward 
Elias 29 with the conditions around several well studied lu- 
minous, high mass protostars. The high gas temperature 
and gas/solid ratios resemble those of relatively evolved 
high mass objects (e.g. GL 2591). However, none of the ice 
band profiles shows evidence for significant thermal pro- 
cessing, and in this respect Elias 29 resembles the least 
evolved luminous protostars, such as NGC 7538 : IRS9. 
Thus we conclude that the heating of the envelope of the 
low mass object Elias 29 is qualitatively different from 
that of high mass protostars. This is possibly related to a 
different density gradient of the envelope or shielding of 
the ices in a circumstellar disk. This result is important 
for our understanding of the evolution of interstellar ices, 
and their relation to cometary ices. 
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1. Introduction 


The general picture of low mass star formation has 
been formed since the 1980 ’s with the availability of in- 
frared and millimeter wavelength broad band photometry 
from the ground, and with the BRAS satellite and KAO ob- 
servatory (e.g., Lada & Wilking 1984; Adams et al. 1987; 
Hillenbrand et al. 1992; Andre et al. 1993). A classifica- 
tion scheme was made, where the continuum emission of 
Class 0 and I objects peaks in the submillimeter and far- 
infrared. These objects are still deeply embedded in their 
accreting envelopes. In the Class II phase, the wind of the 
protostar has cleared its surrounding environment, such 
that it becomes optically visible, and shows H i emission 
lines. The continuum emission of these objects peaks in 
the near-infrared, but there is still significant excess emis- 
sion above the stellar continuum. They are believed to be 
surrounded by optically thick dusty disks. Finally, little 
dust emission remains for Class III objects, when the disk 
is optically thin, and planetary companions may have been 
formed. 

Our knowledge of the physical and chemical state and 
evolution of the material surrounding protostars, has pro- 
gressed with the availability of medium and high res- 
olution spectroscopic instrumentation at near and mid- 
infrared wavelengths (~ 2 — 20 ^m). The progress made, 
is best illustrated by the observations of high mass 
protostars, which are bright and easy to observe. The 
(ro-) vibrational bands of various molecules (CO, H2O, 
CH3OH, silicates) were observed from the ground, re- 
vealing profile variations of the 3.07 /xm, and 4.67 /xm 
H2O and CO ice bands (e.g. Smith et al. 1989; Tielens 
et al. 1991; Chiar et al. 1998). This was interpreted as 
evaporation of the volatile CO ice, and crystallization of 
H2O ice in the molecular envelopes. This heating effect is 
strengthened by the detection of hot gas in 4.6 /xm CO 
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observations (Mitchell et al. 1991). With the launch of 
the Infrared Space Observatory in 1995 (ISO; Kessler et 
al, 1996), it became possible to observe all other molec- 
ular bands in the infrared (Whittet et al. 1996). It was 
shown that high mass protostellar evolution can be traced 
in the gas-to-solid abundance ratios (van Dishoeck et al. 
1996; van Dishoeck & Blake 1998), and the profiles of the 
ice bands, in particular solid CO 2 (Gerakines et al. 1999; 
Boogert et al. 2000). Thus, there is overwhelming evidence 
that thermad processing, i.e. evaporation and crystalliza- 
tion of ices in and around hot molecular cores, plays an 
important role in the evolution of high mass molecular 
envelopes. 

The composition and evolution of the molecular ma- 
terial around low mass protostars are not as well studied. 
It seems unlikely that the molecular material evolves sim- 
ilar to that around high m^LSs protostars. Low mass pro- 
tostars evolve much slower, release less raxiiative energy, 
drive less energetic winds, and form disks. It is not esta^ 
lished whether low mass objects possess hot cores as well, 
and whether the ices survive the process of star formation. 
If (some of) the ices survive, are they included into comets, 
and if so, are the ice structure and composition still the 
same compared to interstellar ices? How important are 
energetic processes, such as cosmic ray bombardment, in 
altering the ice composition on the long time scale of the 
formation of low mass stars? 

To investigate the influence of low mass protostars on 
their molecular envelope, we make an infrared spectro- 
scopic study Elias 29, also called WL 15 and YLW 7 (Elias 
1978, Wilking et al. 1983, Young et al. 1986). On a large 
scale, Elias 29 lies in core E, which is in the south-east 
corner of the 1x2 pc extended compact CO ridge L 1688 
(Loren et al. 1990) in the densest part of the p Ophiuchi 
cloud, at a distance of ~160 pc from the earth (Wilking & 
Lada 1983; Whittet 1974). It is the reddest object found 
in the near-infrared survey of this cloud by Elias (1978), 
without a counterpart at optical wavelengths. For our ob- 
servations, we used Elias’ coordinates (J2000): 

a = 16^27"^09^3 5 = - 24°37'21". 

The overall spectrum of Elias 29 is typical for a heav- 
ily embedded Class I source, probably in a late accretion 
phase (Wilking et al. 1989; Andre & Montmerle 1994; 
Greene & Lada 1996; Saraceno et al. 1996). The embedded 
nature is also revealed by its high extinction, and by the 
cold compact envelope observed at millimeter wavelengths 
(Andre & Montmerle 1994; Motte et al. 1998). Elias 29 
is associated with a molecular outflow (Bontemps et al. 
1996; Sekimoto et al. 1997). With a bolometric luminos- 
ity of ~ 36 Lq (Chen et al. 1995), Elias 29 is the most 
luminous protostar in the p Oph cloud, which makes this 
source very suitable for spectroscopic studies. The rela- 
tively high luminosity, and high bolometric temperature 
(Tboi ~ 410 K) imply an age in the range 0.5-4 10® yr 


(Chen et al. 1995). In the pre-main-sequence evolution- 
ary tracks of Palla & Stabler (1993), this corresponds to 
a star with end mass 3.0-3.5 Mq, Elias 29 might thus be 
a precursor Herbig AeBe star. This classification is how- 
ever uncertain. For example, it has been argued from the 
SED, and the absence of mid-infrared emission features, 
that Elias 29 is a 1 Mq protostar with a large accretion 
luminosity, and a spectral type of K3-4 at the birth line 
(Greene & Lada 2000). 

This Paper is structured as follows. Technical details 
on the ISO infrared observations are given in Sect. 2. All 
the observed emission and absorption features are dis- 
cussed in detail in Sect. 3. Section 3.1 gives a description 
of the continuum shape, and a comparison to other lines 
of sight. The ice composition and thermal history, and the 
silicate band depth with inferred extinction and column 
densities toward Elias 29 are discussed in Sect. 3.2. Then, 
numerous lines of gaseous CO and H 2 O are detected, and 
modeled to derive gas temperatures and column densities 
(Sect. 3.3). The molecular abundances and gas-to-solid ra- 
tios of Elias 29 are compared to a sample of sight-lines, 
ranging from dark cloud cores to evolved protostars. A 
comparison with high mass protostars is made (Sect. 4.1). 
Section 4.2 discusses the origin of the wealth of observed 
emission and absorption features and puts them in a ge- 
ometrical picture, where we review the evidence for an 
extended envelope and an accretion disk. We conclude in 
Sect. 5 with a summary and suggestions for future obser- 
vations. 


2. Observations 


2.1. The 2.S-45 pm spectrum 

A low resolution {R = A/AA = 400), full 2.3-45 pm 
spectrum of Elias 29 was obtained with the ISO ShoH 
Wavelength Spectrometer (ISO-SWS; de Graauw et al. 
1996) during revolution 267 (August 10 1996). The ISO- 
SWS pipeline and calibration files, available in July 1998 
at SRON Groningen, were applied. The spectrum is gen- 
erally of good quality, with well-matching up and down 
scans, and no serious dark current problems, except for 
band 2C (7-12 pm). Here, we found that the up and down 
scans deviate over the silicate band. One scan showed 
good agreement with a ground-based spectrum of Han- 
ner et al. (1995), and we used this to correct the deviating 
scan. Standard after-pipeline steps were applied, such as 
low order flat-fielding, sigma clipping and re-binning (see 
also Boogert et al. 1998). The twelve sub-spectra in the 
2-45 pm range match fairly well at the overlap regions. 
Small correction factors (<15%) were applied to correct 
for the band jumps. 
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At selected wavelength ranges (3-3.6, 4-9, and 19.5- 
28 /im), we also obtained high resolution {R = 1500) 
ISO-SWS grating spectra, in revolution 292 (September 
04 1996). These were reduced similarly to the low reso- 
lution spectrum. We found that the overall shape of the 
spectrum neax 4-5 /xm is quite badly affected by detector 
memory effects, presumably due to the occurrence of scan 
breaks (de Graauw et al. 1996). We corrected for this, 
by applying a wavelength-dependent shift to match the 
low resolution spectrum. This does not affect our conclu- 
sions, since the high resolution spectrum was only used to 
study neirrow features. Also, near 6.9 ^m the scans devi- 
ate significantly because of memory effects. This problem 
is reflected in the large error bars given in this paper, as 
they were derived from the difference between the average 
up and down scans. 

2.2, The 4^-190 fim spectrum 

Elias 29 was observed during Revolution 484 (March 14 
1997) with the ISO Long Wavelength Spectrometer (ISO- 
LWS; Clegg et al. 1996). We obtained 15 scans covering 
the range from 43 /xm to 197 /xm in the low resolution 
mode {R ~200) for a total of 2611 sec of integration 
time. The data were reduced using the Off-Line-Processing 
package (OLP) version 7 and the IS 0-SpectrsJ- Analysis- 
Package (ISAP) version 1.3. The spectra were flux cali- 
brated using Uranus (Swinyard et al. 1996). We And that 
at the ISO-SWS/LWS overlap region near 45 ;xm, LWS 
has 35% higher flux than SWS. This difference is only 
slightly larger than the absolute calibration uncertainties 
of the two instruments, and thus it is doubtful that this 
can be ascribed to the presence of extended emission in 
the larger aperture of ISO— LWS (~ 80^^ versus ~ 25^^). 
We therefore decided to multiply the LWS spectrum down 
with this factor. 

3. Results 


3.1, The spectral energy distribution (SED) 

Elias 29 is only visible at wavelengths l 2 urger than 
~1.5 /xm (Greene h Lada 1996; Elias 1978). Our ISO ob- 
servations show that the continuum emission rises steeply 
between 2-3 /xm, reaches a maximum of XF\—U>xlO~^^ 
W cm“^ at A ^^^5 /xm, and is remarkably flat with 
8 X 10”^® W cm"^ between 20 and 100 pm (Fig. 1). The 
emission has dropped to XFx^ 4 x lO”^^ W cm“^ at 
200 /xm, and by four orders of magnitude at 1300 /xm. Our 
near-infrared spectral continuum fluxes are in excellent 
agreement with broad band fluxes from ground-based ob- 
servations (Elias 1978). Also the ground-based small beam 


10 and 20 /xm observations, as well as the large beam 12 
and 25 /xm IRAS fluxes, match the ISO-SWS observa- 
tion well, thus indicating that at these wavelengths the 
emission is well confined within a region of 8" in diam- 
eter (Fig. 1; Lada & Wilking 1984; Young et al. 1986). 
The reasonable match of the ISO-SWS and LWS spectra 
(Sect. 2.2) indicates that also at 45 /xm the emission is not 
very extended (< 25"). At 100 /xm, however, some large 
scale emission may be present, since the IRAS flux, ob- 
served in a 5.5 times larger aperture, is a factor of 2 larger 
compared to the ISO measurement. 

The observed SED of Elias 29 is different from that 
of massive protostars such as GL 2591, and GL 7009S, 
which peak in the far-infrared (Fig. 2). It has been pro- 
posed that the shape of SEDs is independent of the lu- 
minosity of the central object, and is rather determined 
by the total dust column density along the line of sight 
(Ivezic & Elitzur 1997). In the “standard model” of Ivezic 
& Elitzur, GL 2591 and GL 7009S would have a column 
density corresponding to an Ay/ of several hundred. Elias 
29 must have an Ay < 100, because it does not peak in 
the far infrared. However, the flatness of the SED up to 
100 /xm is not reproduced in these models. A lower column 
density alone thus cannot explain the differences between 
the SED of Elias 29 axid massive protostars. Other factors, 
such as a different density gradient, and the presence of a 
circumstellar disk are probably important. We will discuss 
the structure of Elias 29, in relation to the detected gas 
and ice absorption features, in Sect. 4.2. 

3.2. Ice and dust absorption bands 

Numerous absorption bamds of ices and silicates are 
present in the infrared spectrum of Elias 29 (Fig. 3), such 
that hardly any ’clean’ continuum emission is left. We 
identify each band, derive column densities, and, when 
possible, determine the ice mantle composition and ther- 
mal history. The full spectrum also allows a determination 
of upper limits of abundances for undetected, though as- 
trophysically relevant molecules. The column densities for 
the species discussed below are summarized in Table 1. 

3.2.1. H 2 O ice 

The infrared spectrum of Elias 29 shows ail the vibration 
modes of H 2 O ice in absorption (Fig. 4). We see the O- 
H stretching mode at 3.0 /xm (‘Vi, 1 / 3 ” in spectroscopic 
notation), the 0-H bending mode at 6.0 /xm (“ 1 ^ 2 ”)? 
libration or hindered rotation mode at ~12 /xm (“i^”), 
the combination mode at 4.5 /xm (“3x/l” or “x ^2 -I- i^”), 
and perhaps the lattice mode at ^45 /xm. 

The continuum determination is complicated by the 
large width of all these bands. In accordance with other 
studies (Smith et al. 1989; Schutte et al. 1996; Keane et 
al. 2000), we used single blackbodies to flt the continuum 
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Fig. 1. Spectral energy distribution of 
Elias 29, consisting of ground-based ob- 
servations (A <2.4 ^m; Greene & Lada 
1996), an ISO-SWS spectrum (A =2.4- 
45 /im), and an ISO-LWS spectnim 
(A =45-195 fim). The data point at 
1300 is taken from Andre & Mont- 
merle (1994), which we have connected 
with a dashed straight line to the ISO- 
LWS spectrum, to guide the eye. The 
dotted line is the adopted continuum, 
as determined by blackbody fits and by 
hand. The open circles are ground-based 
and IRAS observations (see text). The 
top inset shows a magnification of the 
10-200 ^m region. 


Fig. 2. Spectral energy distribution of 
Elias 29 compared to the high mass 
protostars GL 7009S (Dartois et al. 
1998a), and GL 2591 (van der Tak et 
al. 1999). The spectrum of the Herbig 
Ae star AB Aur is a compilation of con- 
tinuum observations taken from Man- 
nings (1994), and is further discussed in 
Sect. 4.2. The flux scale of each spec- 
trum has been divided by the values 
given in brackets. 


locally, directly adjacent to the absorption bands. For the 
6 fxm band we took into account that laboratory spectra 
of the bending mode of H 2 O ice show a prominent wing on 
the long wavelength side, extended up to 8 /xm (e.g. Hud- 
gins et al. 1993, Maldoni et al. 1998). We simultaneously 
fitted a blackbody continuum, normalized at 5.1 /xm, and 


a laboratory ice spectrum to the observed flux at 8 fxm 
and the shape of the 6.0 ^m feature (Fig. 5). 

The shape of the 6.0 ;xm H 2 O bending mode is partic- 
ularly sensitive to the ice temperature (e.g. Maldoni et al. 
1998). At higher T in the laboratory, the strength of the 
main 6.0 /xm component decreases at the expense of more 
absorption in the long wavelength wing. In the spectrum 
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Fig. 3. Low resolution {R = 400) mid-infrared ISO-SWS spectrum of Elias 29 with a smooth, global continuum (dotted line), 
rather arbitrarily determined by hand and blackbody fits. The vibrational absorption bands of various molecules are indicated. 


of Elias 29, the wing cannot be seen as a separate feature, 
since at 8 fim it blends with the very deep silicate band. 
However, the observed 6.0 band is relatively sharp, 
and it can only be fitted with H2O ice at T < 80 K, with 
a best fit at T = 40 K (Fig. 5). The excellent fit to the 
6.0 /xm band in Elias 29 indicates that the 5.83 and 6.2 
excess absorptions detected toward several massive proto- 
stars (Schutte et al. 1996, 1998; Keane et al. 2000), are 
not seen in this source (Sect. 3.2.7). 

The observed peak position of the stretching mode 
of H2O ice toward Elias 29 is 3.07±0.01 fim. The short 
wavelength wing is well matched with a laboratory ice at 
T =40 K, as for the bending mode (Fig. 5). The long wave- 
length wing, however, is poorly fitted. It has been realized 
since long that light scattering by large ice grains leads to 
extra extinction on the long wavelength wing (e.g. Leger 

Table 1. Colunm densities of ices toward Elias 29. Non- 
detections are indicated with 3<r upper limits. 


molecule 

N [10*^ cm-2] 

H 2 O 

34(6) 

**C02 

6.7 (0.5) 

^^C02 

0.083 (0.005) 

CO 

1.7 (0.3) 

CH 4 

<0.5 

NH3 

< 3.5 

CH 30 H 

< 1.5 

H 2 CO 

< 0.6 

HCOOH 

< 0.3 

OCS 

< 0.015 

XCN 

< 0.067 


et aJ. 1983). To illustrate this, we calculate the extinc- 
tion cross section for spherical silicate grains coated with 
ice mantles, applying the code given in Bohren & Huff- 
man (1983) and the optical constants of Draine & Lee 
(1984) and Hudgins et al. (1993). Indeed, grains with a 
core-hmantle radius of ~0.6 /xm provide a much better fit 
to the long wavelength wing than small grains do (Fig. 5). 
This effect is unimportant for the 6.0 /xm band since it is 
intrinsically weaker, and the grains are smaller compared 
to the wavelength. In a more realistic approach, a distribu- 
tion of grain sizes, as well as constraints to other observ- 
ables such as continuum extinction, the total grain and 
ice column densities, and polarization need to be tedcen 
into account. Although there is a general consensus that 
large grains need to be invoked (e.g. Leger et al. 1983; 
Pendleton et al. 1990; Smith et al. 1989; Martin & Whit- 
tet 1990), there is no unified grain model yet that obeys all 
the observational constraints (e.g. Smith et al. 1993; Tie- 
lens 1982). Alternative absorbers at the long wavelength 
wing have been proposed, such as H2O.NH3 bondings. As 
illustrated in Fig. 5, this effect is however small at the low 
column density ratio of NH3/H2O < 0.13 toward Elias 29 
(Sect. 3.2.7). A smeJl contribution is also made by absorp- 
tion by hydrocarbons (Sect. 3.2.4). 

The peak optical depth of the 3.0 /xm band is 
1.85i0.08, which is in excellent agreement with the study 
of Tanaka et al. (1990). Using an integrated band strength 
A = 2.0 X lO"^^ cm molecule^S we derive a column den- 
sity of 7V(H20)= 3.0 X 10^® cm"^ for the small grain 
model, and 3.7 x 10^® cm”^ for the large grain model. 
Since at present we can not favor one of these two cases, we 
will assume an average value of iV(H20)= (3.4±0.6)xl0^® 
cm””^ in this paper. The error bar also includes the uncer- 
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X (>^) 

Fig. 4. Optical depth spectrum of Elias 29, assuming the con- 
tinuum indicated in Figs. 1 and 3. The light, thick line is a 
laboratory spectrum of H 2 O ice at T =10 K (Hudgins et al. 
1993). All five H 2 O ice vibration bands can be discerned. 


tainty in band strength, which increases with 10% when 
the ice is heated from 10 to 100 K (Gerakines et al. 1995). 
Note that a column density determination from the 6.0 /zm 
bending mode is more uncertain due to the unreliable con- 
tinuum on the long wavelength side (Fig. 5). At this col- 
umn density of H2O ice, the depth of the other vibrational 
modes is in good agreement with the observed spectrum 
of Elias 29 (Fig. 4). 



X (jitn) 


Fig, 5, Analysis of the 3.0 fim (a) and 6.0 fjtm (b) absorption 
bands of H 2 O ice. In panel a, the thick gray line is a laborar 
tory spectrum of pure H 2 O ice at T =40 K. The dashed line is 
a spectrum of H20:NH3 = 10:1 (T = 50 K). The dotted line 
gives a calculated band profile of a large ice grain (see text). 
Panel b shows the 5-8 /zm spectrum, with laboratory spectra 
at T = 40 K (thick gray line), and at T = 100 K (thin, solid 
line), showing that only low temperatures provide good fits to 
Elias 29. The dashed line shows the assumed blackbody con- 
tinuum. The narrow absorption lines in the observed spectrum 
originate from H 2 O vaj)or (Sect. 3.3). 


3.2.2. CO ice 

The CO ice band at 4.67 /zm in Elias 29 is contami- 
nated by gas phase CO lines from low J levels (Figs. 6 
8ind 10). In particular, the P(l) line lies in the cen- 
ter of the ice band at 4.674 /zm. To study the band 
profile, we subtracted a model for the gaseous lines at 
Tex = 750 K, AT = 5 X 10^® cm“^, and 6 d= 5 km s”^ 
(Sect. 3.3). This increases the ice band width by 0.9 cm“\ 
to FWHM=4.40 cm“^ (0.010 /zm). With a peak position 
of 4.673 /zm (2140.1 cm“^), the CO ice band observed to- 
ward Elias 29 is similar to that of the lunainous protostar 
NGC 7538 : IRS9 (Fig. 6; Tielens et al. 1991; Chiau: et 
al. 1998). The main, narrow component at 4.673 /zm is at- 
tributed to pure solid CO, or CO embedded in an environ- 
ment of apolar molecules. In particular, mixtures with O2, 
at an O2/CO ratio as much as 5 (Elsila et al. 1997; Chiar et 
al. 1998) provide good fits. Mixtures of CO with CO2 are 
generally too broad (Ehrenfreund et al. 1997). While the 
apolar, volatile component dominates the spectrum, both 
Elias 29 and NGC 7538 : IRS9 show evidence for a wing on 
the long wavelength side. This is attributed to CO diluted 
in a mixture of polar molecules such as H2O and CH3OH 
(Chiar et al. 1998; Tielens et al. 1991). Assuming a band 


strength A = 1.1 x 10^^ cm molecule for both the polar 
and apolar components (Gerakines et al. 1995), we derive 
7V(CO ice)=1.7xl0^^ cm”^ with an apolar/polar ratio of 
~8, comparable to NGC 7538 : IRS9. These results are in 
good agreement with the ground-based study of Kerr et 
al. (1993). Although NGC 7538 : IRS9 seems to have a 
larger polar CO component in Fig. 6, this difference may 
merely reflect uncertainties in the continuum subtraction, 
and the fact that the NGC 7538 : IRS9 spectrum is not 
corrected for gas phase CO lines. 

3.2.3. CO2 ice 

The absorption bands of CO2 ice are prominently present 
in the infrared spectrum of Elias 29 (Fig. 3). We see the 
stretching and bending modes at 4.27 and 15.2 /zm respec- 
tively. Not visible in this spectrum is the stretching mode 
of solid ^^C02 at 4.38 /zm, although the high resolution 
spectrum (Fig. 10) shows a hint of its presence. A very 
sensitive observation is presented elsewhere (Boogert et al. 
2000). The ^^C02 bending mode and the ^^C02 stretching 
mode have proven to be very sensitive to ice mantle com- 
position and thermal history. In Elias 29, these bands do 
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X (/xm) 


Fig. 6. The CO ice band on optical depth scale observed to- 
ward Elias 29 (thin solid line). The thick gray line is the spec- 
trum with a gas phase CO model subtracted (see text). The 
dotted line represents the spectrum of the high mass protostair 
NGC 7538 : IRJS9 (divided by 7.5; Tielens et al. 1991) showing 
the similarity of the band profiles. 


not show the narrow substructures seen in many other pro- 
tostars, and attributed to heated polar CO2 ices (Boogert 
et al. 2000; Gerakines et al. 1999). As for the CO ice 
band (Fig. 6), the width and peak position of the ^^C02 
band very much resemble that of the luminous protostar 
NGC 7538 : IRS9. Thus, the CO2 ice toward Elias 29 is 
mixed in with polar molecules, and is not much affected 
by heating. The ^^C02 column density is 22±4% relative 
to H2O ice, which is comparable to the values reported for 
high mass protostars (Gerakines et al. 1999). Finally, we 
derive an isotope ratio of ^^C02/^^C02=81±11 in the ice 
toward Elias 29, which is well within the range found for 
the local ISM (Boogert et al. 2000), 

3.2.4. The 3.47 /xm band 

The long wavelength wing of the deep 3.0 /xm absorption 
band shows a change of slope at 3.38 /xm, indicative of 
a shallow absorption feature (Fig. 7). This feature is £ilso 
detected in an independent ground based study of Elias 29 
(Brooke et al. 1999). For consistency with ground based 
studies, the continuum on each side of the feature was as- 
sumed to start at 3.37 and 3.61 /xm. It must be emphasized 
however, that in particular on the long wavelength side, 
the continuum is poorly defined. Fitting a smooth 6-th 
order polynomial results in an absorption band centered 
on 3.49±0.03 /xm with a peak optical depth of r=0.06 
(Fig. 7). The width is FWHM=120±40 cm’“^ where the 
uncertainty includes the poorly constrained continuum on 
the long wavelength side. Features of similar width and 
peak position have been detected in several massive pro- 



3.2 3.4 3.6 3.8 

X (/xm) 


Fig. 7. a and b. Spectral structure in the long wavelength 
wing of the 3.0 /xm band, a the merged high (R = 1500) and 
low {R = 400) resolution spectra and the assumed polynomial 
continuum (dotted line), b optical depth plot of the detected 
3,47 /xm feature. The gray, thick line represents the ground 
based spectrum of the high mass protostar NGC 7538 : IRS9, 
divided by a factor of 2.3 (Brooke et al, 1999), showing the 
C-H stretch mode of solid CH3OH at 3.54 /xm. This feature is 
absent in the spectrum of Elias 29. 


tostellar objects (Allamandola et al. 1992) and in low mass 
objects and quiescent molecular cloud material (Chiar et 
al, 1996). A likely candidate for this 3.47 /xm band is the 
C-H stretching mode of hydrocarbons. Prom the correla- 
tion of peak optical depths of this feature and the 3.0 /xm 
ice band, it is concluded that the carrier for the 3.47 /xm 
band resides in ices rather than in refractory dust (Brooke 
et al. 1996). We find that with r(3.47 /xm)=0.06 and 
r(3.0 /xm)=1.85, Elias 29 follows this correlation very well. 

3.2.5. CH3OH ice 

In several high mass protostars, the 3.47 /xm band 
is blended with a distinct naurow feature centered on 
3.54 /xm ( Allamandola et al. 1992), This feature is ascribed 
to the C-H stretching mode of solid CH3OH, A direct 
comparison with the high mass protostar NGC 7538 : IRS9 
shows that, although the 3.47 /xm bands have similar 
shapes, the 3.54 /xm feature is absent in Elias 29 (Fig. 7). 
We determine a 3a upper limit to the peak optical depth 
of r(3.54 /xm)<0.036. Scaling with the observed depth and 
column density in NGC 7538 : IRS9 (Brooke et al. 1999), 
then results in an upper limit to the CH3OH ice colunm 
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X (fim) 

Fig. 8. Optical depth plot of the 5-8 fim region of Elias 29, 
after subtraction of an H 2 O ice spectrum at T =40 K as well 
as am H 2 O gas model at Tex = 300 K, N = 2 x 10^® cm"^, 
bu—2.5 km s“^. This figure highlights the 6.85 fim absorption 
feature. The thick gray lines give a compatrison with the mas- 
sive protostars NGC 7538 : IRS9 (top) and S 140 : IRSl (bot- 
tom; shifted down 0.08 along the optical depth axis). 


density 7V(CH30H ice)< 1.5 x 10^^ cm“^, or less than 5% 
of H2O ice toward Elias 29 (Table 3). The other modes 
of CH3OH ice are either much weaker, or aire severely 
blended with the strong H2O and silicate bands (e.g. 
the C-0 stretching mode at 9.7 fim; Schutte et ad. 1991; 
Skinner et aJ. 1992) and thus do not provide better con- 
stradnts on the CH3OH ice column density. Toward other 
low mass objects, and quiescent dark clouds, low upper 
limits have been set to the CH3OH ice abundance as well. 
The CH3OH ice abundance found in massive protostairs is 
generally of the saune magnitude (Chiar et ad. 1996), but 
in a few objects significantly laurger (Dartois et ad. 1999), 
them these upper limits. 

3.2.6. The 6.85 /xm bamd 

Elias 29 is the first low mass protostar in which the 
6.85 /xm absorption band is detected (Fig. 5). After sub- 
traction of the H2O ice band and the gas phase H2O lines 
(Fig. 8), we find that it has a peak optical depth of r ~0.07 
and an integrated opticad depth Tint = 7.8 ± 1.6 cm“^ 
When scaled to the H2O ice column density, the strength 
of the 6.85 /xm band toward Elias 29 is similau: to high mass 
protostars (Keane et ad. 2000). The bamd profile, e.g. the 
shaurp edge at 6.60 /xm, agrees very well with several high 
mass objects, in particular those tracing ’cold’ gas and 
dust (NGC 7538 : IRS9, W 33A, GL 989). It clearly devi- 
ates from warmer lines of sight (e.g. S 140 : IRSl; Fig. 8). 
Thus, in this picture, we find that the material responsible 
for the 6.85 /xm band towaurd Elias 29 is not significantly 


thermally processed. Given the low upper limits to the 
CH3OH ice column density toward Elias 29, only a frac- 
tion of the band, as for high mass objects, cam be explained 
by the C-H bending mode of CH3OH ices (Schutte et al. 
1996). For a detadled bamd profile amadysis amd a discus- 
sion on the origin of the 6.85 /xm band, we refer to Keame 
et al. (2000). 

3.2.7. Upper limits to solid CH4, NH3, H2CO, HCOOH, 
OCS, and ‘XCN’ 

Several solid state species have been detected towaird lu- 
minous protostars, but are absent towaird Elias 29. The 
deformation mode of solid CH4 was detected toward pro- 
tostars, with a peak position at 1303 cm”^ (7.67 /xm), 
and a width FWHM=11 cm"^ (Boogert et ad. 1996; Dau*- 
tois et ad. 1998b). For Elias 29 we cam exclude this bamd 
to a peaik optical depth of r <0.03, corresponding to 
N(CH4)/iV(H20)<1.5%. This 3 <t upper limit is compa- 
rable to the detection in NGC 7538 : IRS9 (Boogert et al. 
1996). 

Solid NH3 was detected by its 9.10 /xm inversion mode 
towaird NGC 7538 : IRS9 (Lacy et al. 1998). Using the 
band strength determined in Kerkhof et ad. (1999), the 
NH3 column density is 13% of H2O ice. Recent detec- 
tions in other highly obscured lines of sight give similair 
(W 33A; Gibb et ad. 2000), or a faetor 2 larger NH3 abun- 
dances (Galactic Center; Chiar et ad. 2000). To find this 
band in the deep silicate feature of Elias 29 we take the 
same approach as Lacy et al., by fitting a local stradght 
line continuum to the wavelength regions 8.52-8.69 and 
9.20-9.55 /xm. As a check, we perform the saune procedure 
to the ISO-SWS spectrum of NGC 7538 : IRS9 (Fig. 9). 
We confirm the detection of Lacy et ad., although the 
peak opticad depth r ~ 0.16 is a factor 2 lower in our 
case. We ascribe this difference to the calibration uncer- 
tadnties of ISO-SWS at this wavelength (Leech 2000). For 
Elias 29, a feature with r ~ 0.06 might be present. How- 
ever, due to the poorly defined long wavelength side of the 
continuum (Fig. 9) amd the ISO-SWS cadibration uncer- 
tadnties, we will assume a conservative upper limit to this 
band of r <0.1. This corresponds to a column density of 
N{NRz)< 3.5 X 101^“^ i.e. iV(NH3)/iV(H20)<13%. 
Other vibrational bands of NH3 do not provide better 
constraints. The equally strong N-H stretching mode at 
2.90 /xm (d’Hendecourt & Allamamdola 1986) is hidden in 
the steep wing of the 3.0 fxxa ice bamd, and there is no sig- 
nificant difference in this region between the laboratory 
spectra of pure H2O ice ’and the mixture H2O:NH3=10:l 
(Fig. 5). A similar problem exists for the N-H deformation 
mode at 6.16 /xm, which is hidden in the long wavelength 
wing of the H2O bending mode (Keane et ad. 2000). A 
feature with a peaik optical depth of r <0.025 would be 
expected here (Sandford & Allamandola 1993). When sub- 
tracting water ice amd vapor absorption, a weadc bamd with 
an optical depth r = 0.03 perhaps remains present at the 
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Fig. 9. a ISO-SWS spectra of the silicate band region of 
Elias 29 and NGC 7538 : IRS9, The dashed line is the local 
continuum for the solid NH3 inversion mode, similar to that 
defined in Lacy et al. (1998). b the residuals after continuum 
subtraction for Elias 29 (solid) and NGC 7538 : IRS9 (dotted). 


expected wavelength (Fig. 8). Given the other positive and 
negative structure in the spectrum, we regard this also as 
an upper limit, however. 

The H20-subtracted 5.0-6. 5 wavelength region 

(Fig. 8) does not show the features detected toward high 
mass protostars (Schutte et al. 1996; 1998; Keane et al. 
2000). At 6.25 ^m (not to confuse with the feature of NH3 
ice at slightly shorter wavelength; see above), a feature 
has been associated with absorption by carbonaceous dust 
(PAH). At 5.83 fim a broad feature has been assigned to 
the C=0 stretching mode of solid HCOOH, and a narrow 
feature of solid H2CO (Keane et al. 2000). Scaling the fea- 
tures observed toward NGC 7538 : IRS9 to the lower H2O 
ice band column density toward Elias 29, one would ex- 
pect peak optical depths T5.83 = 0.06 and T6.25 = 0.03. Our 
spectra indicate upper Umits to these features of r <0.03 
(Fig. 8). Thus, in particular the 5.83 ^m feature toward 
Elias 29 is significantly less pronounced compared to high 
mass protostars. Using the band strengths and typical 
widths given in Keane et al. (2000), we derive 3cr col- 
umn density upper limits of iV(H2CO) < 6 x 10^^ cm”^, 
and AT(HCOOH) < 3x 10^^ cm“^. With abundance upper 
limits of 1-2% with respect to H2O, these aldehydes are 
thus minor ice components. For comparison, toward high 
mass objects it is typically 3% or higher. 

An absorption feature has been detected at 
2042±4 cm“^ (4.90 /zm) in lines of sight toward several 
massive protostars (Palumbo et al. 1997). With a width 
FWHM=23±6 cm”\ it has been ascribed to absorption 
by solid OCS. For Elias 29 this feature is not detected 
with a peak optical depth r <0.01 (3c7), corresponding 
to N{OCS)< 1.5 X 10^5 cm“^ or < 0.05% of H2O ice. 


This upper limit is of the same order of magnitude as the 
detections in W 33A and Mon R2 : IRS2 (Palumbo et al. 
1997). 

Finally, toward several high and low mass protostars a 
feature has been detected at ~2166 cm”^ (4.62 /im) with 
a width FWHM--20 cm'^ (Lacy et al, 1984; Tegler et al. 
1995). This feature is absent in Elias 29, with a peak op- 
tical depth T <0.01 ( 3 < t ). If this feature is caused by the 
C=N stretching mode in ‘XCN’, this corresponds to a col- 
umn density 7V(XCN)< 6.7 x 10^^ cm“^, or less than 0.2% 
of H2O ice (applying A = 3 x 10“^*^ cm molecule”^; Tegler 
et al. 1995), This is considerably less than the detections 
made towaurd high mass objects (e.g. W 33A) and several 
low mass objects (Elias 18; L 1551 : IRS5; Tegler et al. 
1995). This feature has not been detected in the quiescent 
regions of the Taurus molecular cloud (Elias 16; Table 3). 
For a more elaborate discussion on this feature, and the 
proposed carriers, we refer to Pendleton et al. (1999). 

3.2.8. Silicates 

The absorption bands of the Si-O stretching and bending 
modes of silicate dust are prominently present at 9.7 fim. 
and 18 fxm (Fig. 3). We derive a peak absorption opti- 
cal depth of the 9.7 fxm band T9.7 = 1.38 (Fig. 4), which 
is in excellent agreement with the ground-based study of 
Hanner et al. (1995). It is likely that this is a lower limit, 
since the absorption bands have been partly filled in with 
silicate emission from hot dust near the protostar. Model- 
ing of the 9.7 fim silicate band toward Elias 29, including 
emission and absorption, shows that Tg,7 ranges between 
1.51 and 3.38 for optically thick and thin emission respec- 
tively (Hanner et al. 1995). A better fit is obtained for 
optically thick emission. In contrast, for luminous proto- 
stars optically thin emission has been generally assumed. 
Using the relation rg.7 = 1.4 r9.7(obs) + 1.6 (Gillet et al. 
1975; Willner et al, 1982), yields T9.7 = 3.53 for Elias 29. 

For these values of T 9 . 7 , the visual extinction Ay 
ranges between 28 and 65, assuming the standard rela- 
tion Av/t9.7=18,5 (Roche & Aitken 1984). However, these 
limits are likely overestimated (30-50%), because of the 
anomalous extinction curve due to larger grains in the 
p Oph molecular cloud (Bohlin et al. 1978; Martin & Whit- 
tet 1990). Independent extinction determinations, such as 
Ay <48 from the H-K broad band color and Ay <80 
from C^®0 observations (Wilking & Lada 1983), do not 
help to solve this issue. Millimeter continuum observa- 
tions (Andre & Montmerle 1994), and the near-infrared 
J-H color (Greene, priv, comm.), suggest a relatively low 
Ay <30. 

The total hydrogen column density JVh = 
iV(Hl)+2A/'(H2) is closely related to T9.7, and, in 
contrast to the derivation of Ay, the derived is 
not strongly affected by the large grain size in p Oph, 
Applying standard conversion factors for the diffuse 
ISM (Bohlin et al. 1978; Roche & Aitken 1984), we 
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Fig. 10. High resolution (i?=2000) 
spectrum of Elias 29 showing many gas 
phase CO lines^ and the CO ice band 
at 4.67 /xm. The narrow emission line 
at 4.653 /xm is Pf ^ of Hi. Some of the 
narrow structure seen on the long wave- 
length side is due to lines of hot H 2 O va- 
por (see text). For comparison we show 
a low resolution (i^=400) observation, 
shifted along the flux scale for clarity, 
where the gas lines are smeared out over 
the continuum. 


find N}i = 0.5 — 1.2 x 10^^ cm“^, depending on the 
applied rg.y. To be consistent with studies of high mass 
protostars, we will assume in the abundance calculations, 
the value corresponding to optically thin silicate emission, 
i.e. the high limit Na = 1.2 x 10^^ cm”^ (Table 3). 

3,3, Gas phase absorption lines 


is likely a temperature gradient along the line of sight, 
as expected for a protostellar envelope. The LTE assump- 
tion may not apply for the high rotational levels, which 
have high critical densities. Also, the energy levels may be 
pumped by infrared photons, rather than being collision- 
ally excited. Bearing these caveats in mind, we will here 
focus on deriving CO and H 2 O gas column densities and 
temperatures using the LTE models. 


The high resolution 4.00-8.50 /im spectrum of Elias 29 
shows an impressive number of narrow absorption lines of 
gaseous CO and H 2 O (Figs. 5 and 10). We determined 
local continuum points by hand and connected these, us- 
ing a smooth cubic spline interpolation. Then the data 
were converted to optical depth scale, and the absorp- 
tion lines were modeled, using the ro- vibrational spectra of 
gaseous CO and H 2 O described in Helmich (1996). These 
models assume the gas is in Local Thermodyneimic Equi- 
librium (LTE), and has a single excitation temperature 
Tex- The absorption lines have a Voigt profile, and are 
Doppler broadened to a width 6 d (=FWHM/2 vln2). The 
line oscillator strengths are calculated from the HITRAN 
database (Rothman et al. 1992). Finally, the spectrum is 
convolved with a Gaussian to the resolution of our obser- 
vations (i? = 1500 — 2000). Thus, three parameters are 
varied to fit the observed absorption lines: the column 
density N, the Doppler parameter and the excitation 
temperature Tex- Reliable column densities can only be de- 
rived if 6 d is a priori known, which in many studies (like 
ours) is not the case, since the lines are unresolved. At 
low values of the lines become easily optically thick, 
and much larger column densities are needed to fit the 
observed lines, compared to models with high bo values, 
and optically thin lines. 

We emphasize that our assumptions of collisional ex- 
citation, and LTE at a single Tex need not be valid. There 


3.3.1. CO gas 

The 4.4-5.0 pm region shows absorption lines of gas phase 
^^CO, up to rotational quantum number Jiow“33 in the R- 
branch, and Jiow=36 in the P-branch (Fig. 10). The P(l), 
P(2) and R(0) lines are blended with the CO ice hand at 
4.67 pm and the Hi Pf /3 emission line at 4.653 pm. For 
all other absorption lines we determined equivalent widths 
to construct a rotation diagram. A rotation diagram gives 
a first impression of the temperature components present 
along the line of sight, as well as their column densities (or 
lower limits for optically thick lines). For technical details 
on constructing such a diagram we refer to Mitchell et al. 
(1990), and Boogert et al. (1998). The equivalent widths 
were converted to column densities, using the oscillator 
strengths of Goorvitch (1994). For (Fig. 11), we find 
two regimes with very different slopes, corresponding to 
temperatures Trot = 90 ± 45 K and Tot = HOO ± 300 K 
respectively (with 3a errors). However, the slopes of the 
R- and P-branch lines of the hot component are different 
(Fig. 11), resulting in Tot = 1700 ± 420 K when fitting 
to the P-branch lines only. A possible explanation is that 
CO is excited by continuum photons rather than collisions. 
The rising continuum may lead to a higher Tot for the P- 
branch with respect to the R-branch. This effect becomes 
stronger when the photons released after de-excitation of 
R-branch levels are re-absorbed in P-branch levels. Ra- 
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Fig. 11. a rotation diagram of the lines detected toward 
Elias 29 showing the presence of hot and cold gas along the 
line of sight by the different slopes at high and low rotational 
levels, b Xu contour diagram of model fits to the observed ro- 
vibrational spectriun of the R-branch of gaseous CO toward 
Elias 29. x? values are shown for the temperature Tex versus 
CO column density N at constant velocity broadenings 6 d of 
5 km s~^ (dotted) and 10 km s“^ (solid). We only show models 
that provide acceptable fits to the data, i.e. xl < 4- 


diative excitation has also been used to explain the H 2 O 
ro- vibrational spectrum toward Orion BN/KL (Gonzalez- 
Alfonso et al. 1998). The fact that the H 2 O P-branch lines 
are seen in emission for Orion BN/KL, rather than in ab- 
sorption as for CO (and H 2 O; Sect. 3.3.2) toward Elias 29, 
may reflect a different density gradient toward Elias 29, 
such that the photons are not able to escape the enve- 
lope. Additionally, collisionaJ excitation in shocks may be 
of less importance in Elias 29 compared to Orion BN/KL. 
A more careful analysis is needed to discriminate between 
the radiative and collisional excitation models, and alter- 
native explanations, such as non-LTE effects. 




Fig. 12. a Rotation diagram of the ^^CO lines detected to- 
ward Elias 29. Circles are P-branch lines, and triangles are R- 
branch lines. Open symbols refer to ^^CO lines heavily blended 
with ^^CO lines, and are not used to determine the physi- 
cal parameters. The straight line indicates the best fit, with 
a gas temperature Trot = 85 ± 57 K, and a column density 
N(^^CO)=(l.l ± 0.2) X 10^*^ cm“^. b xl contour diagram of 
model fits to the observed ro-vibrational spectrum of gaseous 
^^CO toward Elias 29. xl values are shown for the temperature 
Tex versus CO column density N for constant velocity broad- 
enings 6d=2.5 km s“^ and 6 d= 10 km s”^. Only acceptable fits 
to the data, having xl <3.5, are shown. 


The column densities that we derive from the ab- 
scissa in the rotation diagram are N{CO)^ 1.7 x 10^^ and 
N{CO)= 3.5 X 10^^ cm"^ for the cold and hot CO compo- 
nents respectively. To better constrain the column densi- 
ties and derive more reliable temperatures, one has to take 
into account opticaJ depth effects, using the LTE model 
spectra discussed above. We chose to fit to the frequency 
range 2170-2290 cm“^ ( Jiow > 7 in R-branch), thus mini- 
mizing the contribution from the cold CO component and 
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Fig. 13. The P branch of gas phase CO observed in Elias 29 (top) compared with a well fitting gas model at Tex = 750 K 
(iV = 5 X 10^® cm“^, 6 = 5 km s“^). The bottom panel shows the residual after subtraction of the gas model, which 

contains lines of cold gas. For comparison, a model is plotted (Tex = 50 K, ./V = 2 x 10^^ cm”^, 6 = 2.5 km s“^). 


contamination by lines (see below). We find that 

good fits to these high R-branch lines are obtained only 
for line widths 6 d> 3 km s“^ . Sub-millimeter emission line 
studies indicate 6d=3.6 km s“^ for CO J= 6 — ► 5, but 
much lower values of 6d=1-2 km s“^ for C^®0 J= 1 — ► 0 
smd CS J= 5 4 (Boogert, Hogerheijde, et al., in 

prep.). Indeed, studies of other sources have shown that, 
as a rule, infrared absorption lines are broader than sub- 
millimeter emission lines (van der Tak et al. 1999). Fig- 
ure 11 shows the xi contour diagram of temperature ver- 
sus column density for two values of the line width 6 d= 5, 
and 6d=10 km s”^. The best fitting models have temper- 
atures Tex = 1100i:400 K, in good agreement with the 
rotation diagram. At 6 d= 10 km s”^ the column density is 
well constrained to iV(CO)=(1.3± 0.5) x 10^® cm”^, which 
is a factor of 3 laurger compared to that derived from the 
rotation diagram. Thus at 6 d= 10 km s~^ the lines are 
still somewhat optically thick. At lower 6 d= 5 km s“^, the 
lines become very optically thick, said the column density 
is poorly constrained. Although the best fits with xl< ^ 
have A^(CO)=(8±4)xlO^® cm”^ at Tex = 650±150 K, rea- 
sonable fits sure obtsdned at any N{CO)> 2 x 10^® cm”^ 
for this hot CO gas. 

Sever sd ^®CO lines can be seen in between the ^^CO 
P-branch lines (Fig. 13). At the resolution of our observa- 
tions, the blending with the ^^CO lines hinders analyzing 
the much wesdeer ^®CO lines. But using several well sepa- 
rated lines, we were able to construct a rotation diagram 
(Fig. 12). We find that they result from cold gas at Trot = 
85±57 K (3<7 error), in good agreement with the cold ^^CO 
gas temperature. In the optical thin case, the column den- 


Table 2. Gas phase said H 2 O colunm densities, derived 
with various methods 


Molecule 

Method 

N [10‘« 
cold 

cm 

hot 

‘"CO 

rotation diagram 

0.17 

0.35 

‘"CO 

LTE, 6d= 10 km s”^ 

-- 

1.3±0.5 

‘"CO 

LTE, 6d= 5 km s“^ 

- 

>2 

^®co 

rotation diagram 

9±2‘^ 

- 


LTE, 6d= 2.5 km s“^ 

16±i0“ 

“ 

H 20 

LTE*", 6d= 5 km s“^ 

- 

0.7±0.4 

H 20 

LTE^ 6d=2.5 kms’^ 

- 

2.4±2.1 

H 20 

LTE^, 6d= 5 km s'^ 

< 1 

0.5 

H 20 

LTE^ 6d=2.5 kms’^ 

< 10 

0.5 


^Converted to A^(^^CO) assuming NC‘^C0)/N{^^C0) = 
80 (Boogert et sd. 2000) 

^Single temperature model with Tex =300 K 

^Double temperature model with fixed Nhot = 5 x 10^^ 


sity of this cold component is iV(^®CO)=(l.I ±0.2) x 10^^ 
cm“^. However, the detected ^^CO lines could still be op- 
tically thick. Therefore, we also modeled the ^®CO spec- 
trum, and determine the xl after subtraction of a good 
fitting hot ^^CO gas model (Figs, 13 and 12). In the op- 
tically thick case, such as for 6d=2.5 km s~^, the column 
density can have a wide range A^(^®CO)=(2 ± 1.3) x 10^^ 
cm”*^. Using the isotope abundsoice ratio ^^CO/^®CO=80 
(Boogert et al. 2000), the inferred cold ^^CO column den- 
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Fig. 14. Optical depth plot of gaseous H 2 O lines observed toward Elias 29 (top), compared with model spectra at various 
temperatures. Column densities are AT = 5 x 10^^ cm“^ for all models, and 6d=5.0 km s”^ (Tex = 50 K models have 6d=2.5 
km s“^). The single component T = 300 K and the two component T = 50 + 500 K modeb give equally good fits. The modeb 
with just cold (T = 50 K) or very hot (T = 1000 K) gas clearly do not fit the data. 


sity is thus A^(^^CO)=(16±10) x 10^® cm"^. There is also 
evidence for lines of warm gas ( Jiow > 9), but at low 
significance (< 2a) and no reliable temperature or column 
density could be derived. 

We conclude that the CO gas along the line of sight 
consists of two temperature components, Trot = 90i45 K 
and Trot = 1100 ±300 K. The column density of both com- 
ponents depends highly on the assumed line optical thick- 
ness (Table 2). Until the intrinsic line width is directly 
observed by very high spectral resolution observations, 
we can only give a lower limit of iV(CO-hot)> 2 x 10^® 
cm“^, while iV(CO-cold) is not well constrained, i.e. 
(16 ± 10) X 10^® cm“^. Given that ATh = 1.2 x 10^^ cm“^ 
toward Elias 29, a total gas phase CO column density 
A^(CO)=12 X 10^® cm~^ is expected, assuming that most 
of the gas along the line of sight is molecular and the con- 
version factor N(H2)/N(CO)=5000 applies (Lacy et al. 
1994). Then, the ratio of hot to cold CO gas along the 
line of sight must be at least 0.2. 

3.3.2, H 2 O gas 

We compare the numerous narrow absorption lines de- 
tected in the 5-7.3 /im spectral region of Elias 29 with 
model spectra of H 2 O vapor at various physical conditions 
(Fig. 14). Clearly, the many lines observed at wavelengths 
longer than ~6.55 //m are explained by H 2 O vapor at a 
high temperature (Tex >100 K). On the other hand, the 
relative weakness of the lines observed in the range 6.55- 
6.65 fim imposes a strict upper limit to the temperature 


of this hot gas (T^x <1000 K). To further constrain the 
gas temperature, and the H 2 O vapor column density, we 
determined the xt ^ large number of models. Rea- 
sonable fits to the full 5-7.3 fim range are obtained for 
temperatures of Tex = 350 ± 200 K. The column density is 
constrained to N={7 ±4) x 10^^ cm~^ for low line optical 
depths (6 d> 5). For narrower lines the column density can 
be an order of magnitude larger. 

In a second approach, we test whether both hot and 
cold H 2 O vapor components could be present sJong the 
line of sight, much like the hot and cold CO compo- 
nents. We fitted the regions 5.5-5.8 and 6.55-7.3 /xm, 
which do not contain lines from the lowest rotational 
levels and thus are particularly sensitive to warm H 2 O 
vapor edong the line of sight (Helmich et al. 1996; Dar- 
tois et al. 1998b). The excitation temperature of this gas 
is Tex = 500 ± 300 K, with column densities simil 2 ir to 
that of the single component model. In the high temper- 
ature regime (Tex ^500 K), the modeled line depths in 
the 6. 0-6. 5 /xm region, tracing colder gas, are significantly 
underestimated. To determine the temperature and col- 
umn density of this possible cold component, we fitted 
the sum of a good fitting hot gas model (Tex =500 K, 
yV = 5 X 10^^ cm”^, 6 d= 5.0 km s“^) and a grid of models 
at a wide range of physical conditions to the spectrum of 
Elias 29. Thus, here we assume that the lines of the hot 
and cold gas have different radial velocities and the optical 
depth spectra can simply be added. We find that indeed 
a significant amount of 'cold’ H 2 O vapor, at Tex < 200 K 
may be present (Fig. 14). At Tex < 100 K the column 
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Table 3. Line-of-sight averaged solid and gas phase abundances {N/Nu in units of lO"®) 


species 

Dense Cloud“ 

N7538/9 

W33A 

Elias 29 

GL 2591 

Refs.*" 

H20-ice 

64 

50 

39-143 

28 (8) 

10 

[1], [2], [3, 4], [5], [6] 

-gas^ 

<1 

< 3.1 

<3.6 

>3 

24(3) 

[7],[81,[8].[5],[8] 

t 

6 

0 

17(1) 

8(1) 

3.2 (1.8) 

1.4 ( 0 . 2 ) 

« 1 

[9],[9j,[9],[5],[10] 

-ice [apolar] 

14 (1) 

7 (0.5) 

0.8 (0.1) 

1.2 (0.2) 

« 1 

[9].[91,[9],[5],[10] 

-gas 

< 10 

91 (51) 

143 (32) 

> 67^ 

113 (15) 

[11], [12], [12], [5], [12] 

C02“ice 

12 (3) 

10(1) 

5.2 (0.5) 

5.4 (0.5) 

0.9 (0.1) 

[13],[13],[13],[13],[13] 

-gas 


0.05 (0.01) 

0.08 (0,02) 

<0.06 

0.15 (0.03) 

[8], [8], [8], [8] 

NHa-ice 

- 

7.6 

6.1 

< 3.0 

- 

[14], [3], [5] 

CHaOH-ice 

< 1.8 

2.0 

7 

< 1.3 

4(2) 

[15],[15],[15],[5],[16] 

H2CO-ice 

- 

1.9 

2.5 

< 0.5 

- 

[4], [4], [5] 

HCOOH-ice 

- 

1.1 

0.6 

< 0.3 

- 

[4], [4], [5] 

CH4-ice 


0.8 

0.6 

< 0.4 

- 

[17], [17], [5] 

OCS-ice 

< 0.13 

- 

0.07 

< 0.02 

- 

[18],[18],[5] 

‘XCN’-ice 

< 1.3 

1 (0.3) 

3.6 (1) 

< 0.06 

- 

[19],[19],[19],[5] 

Nh [10*^ cm-*]' 

0.39 

1.6 

2.8 

1.2 

1.7 

[20],[20],[20],[5],[21] 
1 . _ 1 ] 


“ Ice abundances are for Taurus dense cloud toward Elias 16. Gas phase abundances are for p Oph cloud. 
^ References from left to right for each colunm: [1] Chiar et al. 1995; [2] Schutte et al. 1996; [3] Gibb et al. 


2000; [4] Keane et al. 2000; [5] this work; [6] Smith et al. 1989; [7] Liseau & Olofeson 1999; [8] Boonman et 
al. 2000; [9] Chiar et al. 1998; [10] van Dishoeck et al. 1996; [11] Caux et al. 1999; [12] Mitchell et al. 1990; 
[13] Gerakines et al. 1999; [14] Lacy et al. 1998; [15] Chiar et al. 1996; [16] Schutte et al. 1991; [17] Boogert 
et al. 1998; [18] Palumbo et al. 1997; [19] Tegler et al. 1995; [20] Tielens et al. 1991; [21] this work, Fig. 2 
^ All models assume 6d=5 km s~^ for hot gas and 6d=2.5 km s”^ for cold gas 
Total CO ice abundance given as well as the abimdance present in the polar and apolar ice components 
along the line of sight 

® Determined from the 9.7 pm silicate band 


density exceeds the assumed hot H 2 O column density of 
TV = 5 X 10^^ cm”^. For a line width of 6d=5-0 km s"^, we 
find that iV < 1 x 10^® cm“^. At 6 d= 2.5 km s“\ the col- 
umn density of this cold H 2 O gas cannot be constreiined. 

To summarize, the lines in the 5-7.3 ^m range are rea- 
sonably fitted with H 2 O models at Tex ~ 350 ±200 K, and 
iV=(7±4) X 10^^ cm"^ at low line optical depths. For nar- 
rower lines (6 d< 5 km s"^), the column density can be an 
order of magnitude larger. In accordance with the gaseous 
CO along the line of sight, equally good fits are obtedned 
with a two component model, where the cool component 
has Tex < 200 K, and the warmer component Tex >500 K. 
The cool component is then at least as abundant as the 
waurm H 2 O gas. 

4. Discussion 


4 .T Gas and solid state abundances 

We have calculated line of sight averaged gas and solid 
state abundances toward Elias 29, by dividing the col- 
umn densities derived in this paper over the total hydro- 


gen column density Nu = 1.2 x 10^^ cm“^ (Sect. 3.2)^. 
We compare these abundances with a sample of sight- 
lines, spanning the range from dark cloud core to fadrly 
evolved protostars (Table 3). As a tracer of ices in quies- 
cent dark cloud material, we chose the object Elias 16, an 
evolved star by chance located behind the Taurus molec- 
ular cloud (e.g. Whittet et al. 1998). Gas phase CO and 
H 2 O abundances in dense clouds were taken from ISO- 
LWS studies (Caux et al. 1999; Liseau & Olofsson 1999). 
The least evolved protosteu* in our comparison sample is 
NGC 7538 : IRS9. The infrared spectrum of this deeply 
embedded object is characterized by cold ice (Whittet et 
al. 1996), and the gas phase temperatures and abundances 
indicate a very modest hot core (Mitchell et al. 1990; 
Boonman et al. 2000). W 33 A is more embedded than 
NGC 7538 : IRS9, but does have a significant amount of 
warm gas along the line of sight (Mitchell et al. 1990; 
Lahuis & van Dishoeck 2000; Boonman et al. 2000), and 
has a lower abundance of volatile ices (Tielens et al. 1991). 
The most evolved object in our sample is GL 2591. It is 
a typical high mass hot core source, with low ice abun- 
dances and high gas temperatures. All these protostars 
are associated with infrared reflection nebulae, and have 


^ For actual, local, abundances in high mass protostars we 
refer to Boonman et al. 2000 
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well developed high velocity molecular outflows (Mitchell 
et al. 1991; Bontemps et al. 1996; van der Tak et ki, 2000). 
Finally, it is important to note that all the comparison 
protostars are at least three orders of magnitude more lu- 
minous than Elias 29. This allows an investigation of the 
effect of low and high mass star formation on the molecular 
envelopes. An extensive comparison with low luminosity 
embedded objects is at present not possible, because their 
infrared gas and solid state characteristics have not been 
studied in such great detail. 

The H2O and CO ice abundances decrease for the se- 
quence of quiescent dense cloud to NGC 7538 : IRS9, 
W 33A and GL 2591 (Table 3). At the same time, the gas 
phase H2O abundance, the gas phase CO and H2O tem- 
peratures, as well as the gas-to-solid ratios (Table 4), in- 
crease for these objects. All these effects can be explained 
by evaporation of the ice mantles and heating of the hot 
core. It has been suggested that the observed H2O gas 
may also have been newly formed by reactions of atomic 
O and H2 in warm conditions (T > 200 K) in the cen- 
tral hot core or in shocks created by the outflow (e.g., 
van Dishoeck & Blake 1998). However, the total (gas plus 
ice) H2O abundance decreases for the more evolved ob- 
jects, indicating that H2O is destroyed rather than being 
newly formed (van Dishoeck 1998). The low gas phase 
CO2 abundance in all sources indicates that this molecule 
is destroyed even more efficiently after evaporation from 
the grains (Boonman et al. 2000; Charnley Kaufman 
2000). 

In the proposed heating sequence, Elias 29 is placed af- 
ter W 33 A, and before GL 2591. However, the various ice 
band profiles (H2O, CO, CO2, and 6.85 fixn) in Elias 29, 
indicate little thermal processing, resembling very much 
NGC 7538 : IRS9, rather than W 33 A or GL 2591. The 
combination of high gas phase abundances and tempera- 
tures, together with a lack of signatures of thermal pro- 
cessing in the ice bands, as seen in Elias 29, is remarkable 
and is not seen in high mass protostars. Geometric effects 
may play an important role in the evolution of molecular 
envelopes around low mass protostars (see Sect. 4.2). 

Whereas thermal evaporation can explain the abun- 
dance variations of volatiles such as H2O, CO, CO2, NH3, 
and CH4, other mechanisms are needed to explain the vari- 
ations of solid CH3OH, and XCN abundances among the 
sources in our sample (Table 3). It has been widely consid- 
ered that XCN molecules are formed by energetic process- 
ing of icy grain mantles by stellar or cosmic ray induced 
far-ultraviolet radiation, or by bombardment with highly 
energetic particles (e.g. Lacy et al. 1984, Grim Sz Green- 
berg 1987, Allamandola et al. 1988). The high CH3OH 
abundances toward sources with deep XCN bands, and 
the apparent absence of CH3OH towau^d low mass proto- 
stars and dark clouds might suggest that the energetics of 
nearby massive stars is needed to produce CH3OH (Gibb 
et al. 2000). 


Table 4, Gas-to-solid state column density ratios 


Object 

0 

0 

H20« 

C02^ 

rp b 

warm 

Dark Cl. 

<1 

<0.02 

- 

- 

N7538/9 

12 (6) 

<0.05 

0.005 

180 (40) 

W33A 

45 (22) 

<0.11 

0.015 

120 (20) 

Elias 29 

>53 

>0.23 

<0.011 

1000 (500) 

GL 2591 

>400 

2.4 

0.17 

1000 (200) 


^ Determined from Table 3 

^ Temperature of warm CO gas (Mitchell et al. 1990) 


4^2. The structure of Elias 29 

The variety of dust, gas and ice absorption and emis- 
sion components presented here, and in the literature, al- 
lows us to construct an overaJl view of the structure of 
Elias 29. The scale on which the detected hot CO gas is 
present can be constrained when one assumes that the 
pure rotational high-J CO emission lines detected toward 
Elias 29 with ISO-LWS (Ceccarelli et al., in prep.) are 
emitted by the same hot gas. We fit these observed line 
fluxes, by assuming spontaneous, optically thin emission 
from an LTE level distribution, and leaving the size of 
the emitting region as a free parauneter. For the range of 
column densities and temperatures found to fit the CO ab- 
sorption lines (Fig. 11), we find diameters in the range 85- 
225 AU. Thus, the observed hot CO gas may be present in 
a hot core region with the size of a circumstellar disk. The 
gas could be concentrated in a high density photospheric 
layer above the disk. To sufficiently heat it by radiation 
from the central star, the disk needs to flare outwards, 
rather than being flat (e.g. Chiang & Goldreich 1997). We 
cannot exclude however that the gas is present more uni- 
formly in the hot core, at lower densities. It might then 
also be partly heated by shocks from the outflow close to 
the star. A more detailed modeling of the CO emission 
lines, including departures from LTE, optical depth cor- 
rections, and taking into account excitation by radiative 
pumping, are needed to further confine the location of the 
gas phase CO and H2O components (Ceccarelli et al., in 
prep.). 

The lack of signatures of thermal processing in the 
ice bands, locates the ice in a region shielded from the 
central heating source. The ices could be present in a 
foreground cloud, an extended envelope, or a circumstel- 
lar disk seen close to edge-on. Millimeter wave continuum 
observations indicate an extended envelope, concentrated 
on the infrared source (FWHM=17"; 2600 AU; Andre & 
Montmerle 1994; Motte et al. 1998). For a wide range of 
power law model fits to the far-infrared SED, Andre & 
Montmerle find that the envelope mass is 0.1 Mq, with a 
volume averaged dust temperature of typically T =35 K. 
This temperature is high enough to evaporate the most 
volatile, apolar ices, but too low to induce ice crystalliza- 
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tion. Hence, indeed the observed H 2 O, C02» and probably 
'^6.85” fim ices could be associated with this extended en- 
velope. Some of the apolar CO ice has evaporated in the 
envelope after the formation of the low mass protostar, 
as indicated by the significantly lower CO/H 2 O ice ratio 
toward Elias 29, compau-ed to other sight-lines with little 
thermal processing in the ices (NGC 7538 : IRS9, Elias 16; 
Table 3). In this picture, the detected apolar CO ice could 
thus be spatially separate from the other ices, perhaps in 
foreground clouds, or well shielded in a very cold disk. 

Knowledge of the source structure is essential to in- 
terpret the observed solid and gas phase species. For ex- 
ample, if the ice is present in the disk, rather than in the 
envelope, we must see the disk in a neai ^ge-on config- 
uration. Is there independent evidence for the presence 
of a disk surrounding Elias 29 and what would be its 
orientation? The most direct view is provided by lunar 
occultation observations. A central object with diameter 
of ~1 AU emits 90% of the 2.2 /xm continuum emission 
(Simon et al. 1987). The remaining 10% comes primarily 
from an object of 60 AU in diameter, which could be the 
hot part of a disk (T ~ 1000 K). The strongest spectro- 
scopic disk indicator would be the presence of emission 
or absorption of vibrational overtone band heads of CO 
(e.g. Carr 1989; Najitaet al. 1996). The 2.0-2.5 /xm spec- 
trum of Elias 29 does not show these features, in contrast 
to other protostars in p Oph, such as WL 16 (Greene & 
Lada 1996). However, the absence of CO overtone bands 
does not prove the absence of an (inner) disk (Calvet et 
al. 1991). For example, the Herbig Ae object AB Aur does 
not have detected CO overtone bands, while high spatial 
resolution radio continuum and emission line observations 
provide strong evidence for the presence of a circumstellar 
disk around this object (Mannings & Sargent 1997). 

AB Aur is an interesting comparison source, since it 
has the same luminosity as Elias 29 (~40 L©), and the 
SEDs of both objects are remarkably similar (Fig. 2). The 
flatness of the SED in AB Aur is well reproduced in flar- 
ing disk models, where the dust in the outer parts of the 
disk is more efficiently heated than in flat disks (Chiang & 
Goldreich 1997). The disk is optically thick up to 100 /xm, 
and becomes opticadly thin at longer wavelengths where 
the SED drops steeply (e.g. van den Ancker et al. 2000). 
The similarity of the SEDs does however not necessarily 
imply that Elias 29 is dominated by an optically thick disk 
as well. A flat SED could also be produced by the enve- 
lope, if it has a shallow power law density profile (index 
~ 0.5; Andre & Montmerle 1994). This density profile is 
remaurkably flat compared to high mass protostars (van der 
Tak et al. 2000; Dartois et al. 1998a), and other low mass 
protostars (e.g. Hogerheijde & Sandell 2000). Finally, flat 
energy distributions are also created by the combination of 
a disk and envelope. Here, the heated envelope irradiates 
the outer parts of the disk (Natta 1993). 

Without direct high resolution imaging, it is difficult 
to discriminate between these models. Assuming a given 


model however, the present observations put some con- 
straints. In the disk scenario, its orientation would have 
to be closer to edge-on than face-on to explain the ab- 
sorption line spectrum of Elias 29 (Chiang & Goldreich 
1999). In these models, an inclination larger than ^70*^ 
can however be excluded, because this would give an SED 
that peaks in the far-infrared, in contrast to what is ob- 
served for Elias 29. Also, if the disk were edge-on, a higher 
absorbing column, perhaps an order of magnitude larger 
than the observed Nu ~ 1.2 x 10^^ cm”^ (Sect. 3.2) would 
be expected (Sekimoto et al, 1997). An independent mea- 
sure for N}i and the disk inclination is provided by the 
hard X-ray flux and spectrum, arising from hot gas in the 
magnetosphere. For Elias 29, a high Nn ~ 2 x 10^^ cm“^ is 
observed during X-ray flares, but Nn is a factor of 5 lower 
in quiescent phases (Kamata et al. 1997). Perhaps the X- 
ray flares are formed low in the magnetosphere, and in the 
relatively high inclination of the disk, they trace higher 
column densities compared to X-rays formed in quiescent 
phases higher in the magnetosphere. 

5. Conclusions and future work 

The 1.2-195 /xm spectrum of the low mass protostellar 
object Elias 29 in the p Ophiuchi molecular cloud shows a 
wealth of absorption lines of gas and solid state molecules. 
Hot CO and H 2 O gas are detected (Tex >300 K) at rather 
high abundances, on scales of not more than a few hun- 
dred AU. The ice abundances are relatively low. In this 
respect, Elias 29 resembles luminous protostars with sig- 
nificantly heated cores, such as GL 2591. However, none of 
the many ice bands that are detected, i.e. firom H 2 O, CO, 
CO 2 , and the 6.85 /xm band, shows outspoken signs of 
thermal processing. Again in comparison with luminous 
protostars, Elias 29 now resembles less evolved objects, 
such as NGC 7538 : IRS9. Our combined gas and solid 
state analysis thus shows that high and low mass pro- 
tostars heat their molecular envelopes in different ways. 
This may be related to their different structure, such as 
the presence of a circumstellar disk in low mass protostars. 
The hot gas of Elias 29 could be present on the surface 
of a flaring disk, which is efficiently heated by the central 
star. The ices toward Elias 29 must be well shielded in a 
circumstellar disk seen close to edge-on, or far away in the 
envelope. 

Does this imply that in general the ices in the disks 
or outer envelopes of low mass protostars remain unal- 
tered, both in composition and structure, during the pro- 
cess of star formation? Are these ices the building blocks 
of the early solar system and are they preserved in present 
day observed cometaury nuclei? To date, no Class I proto- 
star has been found with strong signs of crystalline ices 
(Boogert et al. 2000). On the other hand, the presence of 
crystalline ices and silicates has been reported in several 
isolated, less embedded Herbig AeBe objects (Malfait et 
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al. 1999). This research needs to be extended to a larger 
sample of low mass protostars, in a range of evolutionary 
stages and luminosities. Furthermore, it is essential for the 
interpretation of the gas and solid state characteristics to- 
ward Elias 29 that the presence of a circumstellar disk, 
and its inclination are determined by future high spatial 
resolution infrared or millimeter continuum observations. 
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